Interaction of vector solitons with a nonlinear interface 
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We develop the analytical method of field momenta for analyzing the dynamics of optical vector 
solitons in photorefractive nonlinear media. First, we derive the effective evolution equations for 
the parameters of multi-component solitons composed of incoherently coupled beams and investi- 
gate the soliton internal oscillations associated with the relative motion of the soliton components. 
Then, we apply this method for analyzing the vector soliton scattering by a nonlinear interface. 
In particular, we show that a vector soliton can be reflected, transmitted, captured, or split into 
separate components, depending on the initial energy of its internal degree of freedom. The results 
are verified by direct numerical simulations of spatial optical solitons in photorefractive nonlinear 
media. 
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I. INTRODUCTION 

Spatial solitons are self-trapped optical beams for 
which diffraction divergence is exactly compensated 
by nonlinear self- focusing Q. Such solitons possess 
many important properties which are believed to find a 
wide range of potential applications in optical scanning, 
switching, and processing devices. Spatial solitons can 
be generated in photorefractive crystals at low powers 
(of order of fiW) and, therefore, such materials look 
promising from a viewpoint of their use for the soliton 
generation in all-optical information circuits. 

One of the important properties of spatial solitons in 
photorefractive media is their stability against collapse 
and break-up even in a bulk medium due to the satu- 
rating nature of nonlinearity. Spatial solitons in (1+1)- 
dimensional systems can be realized for the light prop- 
agation in weakly guiding planar structures, where self- 
focusing properties of the core material lead to the light 
localization along the plane of a linear guiding structure. 

Usually, a spatial soliton is viewed as a lowest-order 
(fundamental) mode of an effective waveguide it excites 
in a nonlinear medium. Moreover, higher-order modes 
of the induced waveguide, that are associated with an 
effective soliton-induced potential possessing more than 
one energy level, can be excited and such modes can 
be treated as the soliton internal modes. Such inter- 
nal modes are responsible for oscillations of the soliton 
amplitude and width, and they have already been stud- 
ied in a number of theoretical papers H ||. A much 
more interesting situation takes place in the case of the 
so-called vector spatial solitons, i.e. multi-component op- 
tical solitons which consist of several mutually incoherent 



(or orthogonally polarized) co-propagating beams. This 
problem got a little of attention in the past ||, but it 
attracts much attention in view of recent experimental 
studies of multicomponcnt optical solitons. 

The aim of the present paper is twofold. First, we 
study the internal dynamics of (l+l)-dimensional vector 
solitons, associated with the transverse relative oscilla- 
tion of the centers of mass of the partial soliton compo- 
nents. Second, for the first time to our knowledge, we 
study the scattering of a vector soliton by a nonlinear in- 
terface and demonstrate that the scattering process can 
be modified dramatically by initial excitation of the soli- 
ton internal modes. 

As an example, we consider the dynamics of solitons 
in photorefractive nonlinear media taking into account 
saturation of the medium nonlinearity. A standard way 
for studying the excited states of solitons is the anal- 
ysis of the corresponding linear eigenvalue problem for 
small perturbations. However, in the case of nonintc- 
grable systems, such excited eigenstates can be found 
only by rather complicated computer simulations (see, 
for example, Ref. |(| and references therein). In contrast, 
here we suggest the method of momenta and applied it 
to analize the scattering problem. As a matter of fact, 
such a method corresponds to the so-called quasi-classical 
approximation in quantum mechanics. Regardless of the 
fact that the method of momenta leads to the continuous 
spectrum of the eigenstates, it allows us to consider large 
perturbations (including the vector soliton splitting into 
individual components), since this approach is dealing 
with average characteristics of a vector soliton. 



II. VECTOR SPATIAL SOLITONS. 
THE METHOD OF MOMENTA 
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We consider spatial vector solitons consisting of N in- 
coherently coupled components, which are described by 
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the following set of Schrodinger-like coupled nonlinear 
equations (see. e.g. Ref. [Q and references therein): 



dx 



dz 2 



ejvi (I) = 0, 



(1) 



where x and z are the propogation and transverse co- 
ordinates, respectively, (j = 1, ...JV) are the dimen- 
sionless slowly varying amplitudes of the beam compo- 
nents, Xj are the partial propagation constants, £nl (I) 
is the nonlinear correction to the dielectric permittiv- 
ity, / = Y]j—i Ij is the total intensity of the beam, and 

Ij = \^j\ 2 are the partial intensities of the soliton com- 
ponents. We assume that the nonlinear correction to the 
dielectric permittivity depends only on the total inten- 
sity of the light, /, which can be written approximately 
as follows [pi: 



e NL {I) = I/{1 + si) 



(2) 



where s is the saturation parameter. A ground state 
of a vector spatial soliton corresponds to the station- 
ary solution of eq. ([!]) with d/dx = in which tyj = 



(z -Zj), where ^ (z) are real localized functions 
and Zj is the center of mass of the j-th component. As 
a rule, the values Zj for the vector solitons in the ground 
state are equal to each other and coincide with its center 
of mass. The dynamical behaviour of an excited vector 
soliton, such as oscillation of its amplitude and width, 
oscillation of the centers of mass of its components Zj, 
or superposition of these oscillations, depends on initially 
excited intrinsic degrees of freedom. Here, it is important 
to note an analogy between a vector spatial soliton and 
a bound state of quantum particles. Indeed, the differ- 
ence lies only in self-consistent nonlinear nature of the 
potential that provides the confinement for soliton com- 
ponents. 

To study excited states of vector solitons we use the 
method of momenta that has been initially developed for 
scalar solitons of a Kerr- type nonlinear medium |9| . Here 
we generalize this aproach and technique to study multi- 
component spatial solitons in photorefractive media. 

The method of momenta describes the electromagnetic 
beams in terms of integral field characteristics, or field 
momenta 101. The first momentum is defined as 



-,(<>) 



I j dz. 



(3) 



It has the obvious physical meaning of the energy flux in 
the j-th soliton component. The second momentum is 
defined as 



zlj dz. 



(4) 



Here the value Zj = Zj/Pj is the center of mass of the 
j-th beam component. 

The purpose of this section is to derive a close set of 
effective ordinary evolution equations for the momenta 



(Q, H) . Similar procedure is used in semiclassical descrip- 
tion of quantum systems (see, e.g. Ref. 0). Let us 
multiply Eq. (|l|) by ^* to obtain an equation for Pj , and 
by zSfrj to obtain an equation for Zj, and then subtract 
the corresponding complex conjugated equations. The 
integration of the resulting equations over z yields 



dPj_ 
dx 



0. 



dZj 
dx 



\&* 1 

1 dz 



* 3 



dz. 



(5) 



(6) 



Equation (^) means that the energy flux in each indi- 
vidual soliton component is a conserved quantity. Thus, 
there is no energy exchange between the components. 
After differentiation of Eq. (^|) by x and subsequent in- 
tegration over z, we obtain 
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dx 2 



de NL 
dz 



dz. 



(7) 



The physical meaning of the Eqs. (Q) comes from an 
analogy with the geometrical optics: the Eqs. (]?]) deter- 
mine optical rays that correspond to the partial beams 
forming a vector soliton. From the set of equations (Q) 
one can see that the center of mass of the whole beam 
moves along a straight line: 



d 2 Z c 
dx 2 



dx 2 



Z\P\ + ZiPi 



+ Z N P N 



Pl+P-2 



P 



N 



0. (8) 



The set of Eqs. (Q) is not closed. It cannot be used di- 
rectly since, in order to trace the dynamics of the soliton 
components, we should know the electric field at each 
point (x, z) in space, i.e. effectively complete solutions 
of the Eqs. ([!]) are required. In this sense Eqs. (Q) are 
equivalent to Eqs. ([!]), and so far we have not gained 
anything in addition. There are two ways to derive a 
closed set of equations from Eqs. (0). The first one 
is to obtain equations for higher order momenta which, 
nonetheless, does not guarantee that the system of equa- 
tions will become closed. A simpler way, which allows us 
to obtain the desired result, is to employ a trial function 
to be substituted into Eqs. (f7j) instead of the unknown 
functions tyj. It is clear that, for small deviations of 
the soliton constituents from the lowest-order (ground) 
equilibrium state, we can choose the trial functions in 
the form of unperturbed stationary soliton components. 
Moreover, since within Eqs. (]?]) we are operating with 
the average (integral) characteristics of the soliton, one 
can guess that the result of the integration will not be 
sensitive to variations of the transverse soliton structure, 
even for large deviations from the ground state. This is a 
rather strong statement, and it will be verified below by 
means of direct numerical simulations of the Eqs. (m). 
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FIG. 1: (a) Effective potential of the internal motion of the 
soliton components, (b) transverse structure of the soliton 
with identical components, (c) solid - trajectories of the cen- 
ters of mass of the beams calculated by the method of mo- 
menta; dashed - the same values found by numerical simula- 
tions of the Eqs. ([j]). 



III. TWO-COMPONENT SOLITONS 

To demonstrate the internal dynamics of the soliton 
excited states, we consider a particular case of two- 
component vector solitons. As discussed above, it is nat- 
ural to choose the trial function for the Eqs. (0) in the 
form of a stationary solution ty^\z — Zj). According to 
Eqs. ([7]), the relative shift between the soliton compo- 
nents, A = Z2 — Zi , is governed by the following equation: 



(9) 



where P = P1P2/CP1 + P2) is the reduced mass, and 

de NL [h(z) + h{z + A)] 



/OO 
I 2 (z + A)- 
-00 



dz 



dz . 
(10) 



Let us now introduce the effective potential of the beams 
interaction as follows 



F(A) = 



dV(A) 

d~A ' 



(11) 



We have now come to the simple mechanical model 
for the internal dynamics of the two-component soli- 
ton, which describes two interacting quasi-particles in a 
bound state. This state can be thought of as a "photon 
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FIG. 2: (a), (c) Structure of the soliton components at s = 
0.3, Ai = 1, Aa = 0.88 and Ai = 1, Aa = 0.9, respectively, (b), 
(d) Effective potentials of the internal motion for (a) and (c) 
solitons, respectively. 



molecule". Fig. [I] shows the effective potential V(A) for 
numerically calculated one-hump soliton with two identi- 
cal components. The minimum of the effective potential 
corresponds to the stationary soliton in a ground state. 
The internal oscillations of excited soliton components 
deducted from the solution of Eq. (||) (shown in Fig. |l| 
(c) by solid lines) are compared to the numerical solu- 
tion of cqs.(0) (dashed lines in Fig. [j] (c)). The damping 
of the oscillations occurs due to radiation losses of the 
energy from the soliton originating from bends of the 
soliton component trajectories, which are not taken into 
account in the method of momenta. 

The next example deals with the two-hump solitons 
with the transverse structures of the components shown 
in Fig. H (a,c) |Tjj| . One can see that, when the ampli- 
tude of the two-hump component is significantly smaller 
than the amplitude of the one-hump one, the stationary 
state with coaxial constituents is unstable because of the 
maximum of the potential at A = (see Fig. || (b)). It 
is interesting to note that, in this case, the ground equi- 
librium state of the soliton corresponds to the soliton 
components with displaced centers. Such (1+1)D asym- 
metric vector solitons are well known from the theory of 
the two-component Kerr solitons described by Manakov 
model (s — > in eqn. |J) |l2| and for the partially coher- 
ent multi-component solitons in photorefractive medium 
p3[ . However, when the amplitudes of the components 
are comparable, the coaxial state becomes stable again 
(see Fig. ^ (d)). The dynamics of the unstable coax- 
ial soliton state calculated by solving eqs.Q, as well as 
cqs.(Q), is shown in Fig. ||. One can see that, in this 
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mass of the components we have: 
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FIG. 3: Dynamics of a two- hump soliton [shown in Fig. |^ 
(a)] dynamics: solid - trajectories calculated by the method 
of momenta; dashed - the same values found by numerical 
simulations of the Eqs. (jjj). 



case, the agreement between the method of momenta 
and numerical simulations is much worse than for the 
one-hump solitons because of the strong reshaping of the 
two-hump beams. This reshaping, along with the radia- 
tion processes, violates the assumptions of the analytical 
method. 



IV. INTERACTION OF VECTOR SOLITONS 
WITH A NONLINEAR INTERFACE 

The method of momenta looks promising for treating 
the problem of vector soliton dynamics in inhomogeneous 
media. In this section we consider interaction of a two- 
component vector soliton with an interface between two 
photorefractive media. Both linear and nonlinear prop- 
erties of the materials are discontinuous on the interface. 
The problem of one-component photorefractive soliton 
interaction with nonlinear interface had recently been 
studied in |Q, soliton interaction with interface in the 
Kcrr-likc medium can be found in p), |l§]. Let us rep- 
resent the coordinate dependences of the partial prop- 
agation constants and nonlinear correction to dielectric 
permittivity by the following form: 



Al,2(» = Ai, 2 + S\i >2 ■ l(z), 
ejvi(^) = £nl(I) ■ [1 + 5e NL ■ l(z)] 



(12) 



where l(z) 



is the unit step function. By 



l,z > 
0,z < 

the method of momenta one can obtain that the energy 
fluxes in both beams are conserved. For the centers of 



[2 + Se NL -l(z)}Ii,2 



de NL (T) 



~d^ ~ J_ x L " ^ UCiVi ' ^ n '^~dz~ dZ ~ 

-5t NL Il,2{-Zi t2 ) + 5Al,2/l,2(— ^1,2) x 

xcjvl [hpi-Ztp) + h,i (-^2,1)] 

(13) 

All the terms in the Eqs. ([l3|) have clear physical mean- 
ing. The terms that do not contain the jumps of the 
medium parameters represent the internal forces acting 
between soliton components. They are responsible for 
an internal dynamics of the soliton and, in the homo- 
geneous case, can be reduced to the equation (^|). The 
terms proportional to <5Ai.2 and Ssnl are the external 
forces perturbing the internal oscillations of the soliton 
components at the collisions of quasi-particles with an 
interface. 

To derive the closed set of equations describing the in- 
teraction of a vector soliton with an interface we substi- 
tute into Eqs. ( |l3| ) trial functions in the form of the sta- 
tionary unperturbed soliton component (see Section II). 
One can expect that this ansatz will be sufficiently accu- 
rate whenever transverse structures of the beams are not 
strongly affected by large discontinuities of the medium 
parameters. We will verify this assumption by compar- 
ing direct computer simulation of the wave equations (Q) 
(in the case of inhomogeneous media) with the results of 
quasi-particle theory obtained with the method of mo- 
menta. In what follows, we study the following effects 
of soliton interaction with an interface: (a) excitation of 
internal degrees of freedom; (b) the change of the type 
of soliton interaction with an interface (i.e. reflection, 
transmission and capture) as function of the initial inter- 
nal energy of the vector soliton launched under different 
initial conditions; (c) the splitting of pre-excited soliton 
into individual components at the collision with an in- 
terface. It should be noted, with regard to the point (a) 
above, that a spatial soliton consisting of identical com- 
ponents, being initially in ground state, does not become 
excited at the collision with an interface. This happens 
because, due to identical character of interaction of ev- 
ery component with interface, the coaxial structure of 
the soliton state is preserved. 

The results that we present in this section are obtained 
for the one-hump soliton shown in the inset in Fig. |^ (a) 
and interface parameters arc as follows: i5Ai,2 = —0.15; 
Scnl = 0.2; ds — 0. The solutions of the ordinary dif- 
ferential equations ( ^3|) are shown in Fig. |]. Excitation 
of internal oscillations after soliton reflection from an in- 
terface is demonstrated in Fig. [| (a). Soliton dynamics 
near the interface depends dramatically on the initial en- 
ergy level of the internal oscillations of the components. 
In Fig. U (b) soliton components were launched to the 
interface at the same angle (dZi^/dx — —0.5) and from 
the same position of the common center of mass. The 
only adjustable parameter, which determines the initial 
internal energy of oscillations, was the relative shift be- 
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FIG. 4: (a) Excitation of the internal oscillations of the 
soliton at the reflection from the interface, inset shows the 
transverse structure of the soliton, (b) soliton interaction with 
the interface at different values of A: 1 - A = 0.0, 2 - A = 0.4, 
3 - A = 0.8, 4 - A = 1.2, 5 - A = 1.4, (c) A = 2.2, beams 
make up a bound state without interface, (d) beams are split 
at the interface. 
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FIG. 5: Numerical simulations: (a) excitation of the internal 
oscillations of the soliton at the reflection from the interface, 
(b) soliton interaction with the interface at different values of 
initial value of A, (c) separation of the beams at the interface, 
A = 2.2. 




FIG. 6: Numerical simulations: strong reshaping of the two- 
hump soliton at the interaction with an interface takes place, 
surface and contour plots. Initial conditions are taken in the 
form of the soliton shown in Fig^ ^ (a). Parameters of the 
interface are the same as in Fig. 



tween components A. One can see that, depending on 
the value of A, the transmission (curve 5 - A = 0.12), re- 
flection (curves 2,3,4 - A = 0.4, 0.8, 1.2, respectively) and 
capture into the unstable nonlinear surface mode (curve 
5 - A = 0.0) take place. 

It is interesting to note that the character of soli- 
ton interaction depends upon the parameter A non- 
monotonically, which can be explained by considering the 
phase of the internal oscillations at the interface. This 
phase depends on the initial value of A by virtue of an- 
harmonicity of internal oscillations in the effective poten- 
tial well. Fig. H demonstrates the results of numerical 
calculations of equations ([!]), in which Ai^ and ejvz, are 
taken according to eq. (12). Excitation of the internal 
degree of freedom takes place at the reflection from the 
interface Fig. || (a). The calculations shown in Fig. || 
(b) were carried out for the same initial conditions as for 
curves 1,2, and 3 in Fig. |J (b). The higher values of the 
initial relative shift A between individual components do 
not result in the same behaviour of the trajectories as 
m Fig. |(b). It can be easily explained by the strong 
damping of oscillations due to the radiation effects, which 
takes place for large values A. As a result, the amplitude 
of the internal oscillations is small at the point where soli- 
ton meets the interface and the outcome of interaction is 
very similar to the one for the smaller initial values of 
A. To conclude, one can check that method of momenta 
leads to the satisfactory agreement with the direct nu- 
merical solution of nonlinear Schrodinger equations for 
the values of A up to the half width of the soliton. 

The splitting of the vector soliton into individual com- 
ponents is shown in Fig. |[ We have revealed that the 
splitting occurs at the angles of incidence of the soliton 
beam close to the angle of total internal reflection and for 
large enough values of initial relative shift between com- 
ponents (A = 2.2, dZi/dx = dZ 2 /dx = -0.5). Without 
interface, the soliton components make up a bound state 
(Fig. H (c)). At the presence of interface the method of 
momenta predicts that splitting into components takes 
place (Fig. |(d)). This effect can be treated in terms of 
the quasi-particle theory. Indeed, when the depth of ef- 
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fective potential well is of the same order that the energy 
level of the bound state, than even a small perturbation 
(at the interface) of this state can lead to its destruction. 

We have obtained the similar splitting effect [shown in 
Fig. H (c)] by solving the equations (Q) numerically for 
the initial conditions used in the method of momenta. 
In numerical simulations, we can observe the separation 
of the beam components, but in the case of such large 
values of the initial transverse shift, the value of the first 
momentum (|J) fails to determine the position of the cen- 
ter of mass of the localized beams. Significant part of 
the energy is transferred into radiation, and the integra- 
tion in (^) gives us the coordinate of the common center 
of mass of both localized and nonlocalized waves. It is 
clear that the less the radiation losses are, the better the 
position of the localized beam is described by (|j) . 

Finally, it is neccessary to note that the dramatic re- 
shaping of the multi-hump solitons at the interaction 
with the interface takes place. Numerical simulations 
show that these complex beams can be transformed into 
a set of single hump solitons and a considerable part of 
the electromagnetic energy goes into the radiation even 
for a very small change of the parameters of the medium 
at the interface (see Fig. [}]). This reshaping makes it 
impossible to apply the method of momenta to the prob- 
lems of multi-hump soliton dynamics in the inhomoge- 
neous media. 



V. CONCLUSIONS 

We have developed the method of momenta to ana- 
lyze the dynamics of vector solitons consisting of two (or 
more) incoherently coupled components, in both homoge- 
neous and inhomogeneous photorefractive nonlinear me- 
dia. First, we have demonstrated that this method pro- 
vides an effective analytical tool for the study of internal 
oscillations of vector solitons. Second, we have studied 
the scattering of a vector soliton by a nonlinear interface, 
and we have demonstrated that this type of the soliton 
interaction can lead to the excitation of the soliton's in- 
ternal oscillations. In particular, we have demonstrated 
that the same vector soliton incident onto an interface 
can become reflected, transmitted, captured by the inter- 
face, or even split into its constituents, depending solely 
on the energy of the initially excited internal oscillations. 
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